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OTIC  QUALITY  INSPECTED  3 


PREFACE 


This  report  documents  work  done  to  complete  USAFETAC  Project  #910108,  which  was  in 
response  to  a  support  assistance  request  from  the  HQ  Air  Weather  Service  Product  Improvement 
Division  (AWS/XTX).  USAFETAC  was  asked  to  estimate  spatial  and  temporal  errors  inherent 
in  interpolating  weather  observations.  Error  estimates  were  to  be  stratified  by  time  of  day, 
season,  and  geographical  region.  The  results  were  to  be  used  in  assessing  limitations  on  the 
accuracy  of  weather  observations. 

Nearly  all  weather  support  products  depend  on  the  fundamental  accuracy  of  the  basic  weather 
observation.  Observational  accuracy  depends,  at  least  in  part,  on  frequency  and  spacing.  In 
many  parts  of  the  world,  surface  observations  are  only  taken  at  3-hour  intervals.  Weather  station 
spacing  is  irregular,  and  there  may  be  less  than  one  reporting  station  every  200  km. 

When  a  weather  observation  is  needed  at  a  specific  time  and  place  (as  in  climatological  analysis), 
it  is  usually  necessary  to  interpolate  by  using  the  data  closest  to  that  time  from  the  nearest 
available  weather  station.  There  are  errors  inherent  to  this  procedure,  and  users  of  such  data 
normally  require  estimates  of  those  errors.  This  study  provides  these  error  estimates  as  functions 
of  time  and  distance  for  surface  weather  observations  of  temperature,  wind  speed,  and  wind 
direction.  After  the  methodology  is  explained,  results  are  provided  in  appendices  as  probability 
data  in  a  series  of  charts  that  show  percentiles  of  cumulative  distribution  of  changes  in  tempera¬ 
ture,  wind  speed,  and  wind  direction  as  functions  of  distance  or  time. 

The  authors  wish  to  acknowledge  Lt  Col  Roger  C.  Whiton  (USAFETAC/DN)  and  Mr  Charles  R. 
Coffin  (USAFETAC/DNY)  for  their  contributions  to  this  work.  A  special  acknowledgment  is 
extended  to  Mr  Albert  R.  Boehm  of  Hughes  STX  Corporation  (under  contract  to  the  Phillips 
Laboratory’s  Geophysics  Directorate)  for  recommending  the  approach  used  in  this  study  and  for 
providing  several  helpful  suggestions. 
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1.  INTRODUCTION 


1.1  Purpose.  Nearly  all  weather  support  depends  on  the  fundamental  accuracy  of  the  weather  observation. 
Observational  accuracy  depends,  at  least  in  part,  on  frequency  and  spacing.  In  many  parts  of  the  world,  surface 
observations  are  only  taken  at  3-hour  intervals.  Weather  station  spacing  is  irregular  and  there  may  be  less 
than  one  reporting  station  every  200  km.  Therefore,  when  a  weather  observation  is  needed  at  a  specific  time 
and  place  (as  in  climatological  analysis),  it  is  usually  necessary  to  use  data  closest  to  that  time  from  the  nearest 
available  weather  station.  There  are  errors  inherent  in  this  approach,  and  users  of  the  observational  database 
usually  require  estimates  of  those  errors.  The  purpose  of  this  study  is  to  provide  such  estimates  as  functions 
of  time  and  distance  for  surface  weather  observations  of  temperature,  wind  speed,  and  wind  direction. 

1 .2  Approach.  The  spatial  and  temporal  correlation  of  surface  observations  is  a  function  of  many  variables, 
including  distance,  direction,  climate,  geography,  season,  time  of  day.  Although  there  are  no  limits  on  the 
types  of  stratification  possible,  this  study  focuses  primarily  on  distance  and  its  temporal  counterpart, 
observation  age.  The  probability  that  a  temperature  at  a  distances  from  a  given  location  (or  a  time  t  from  a 
given  time)  falls  within  a  specified  range  will  be  computed.  These  probabilities  will  be  stratified  by 
climatological  region,  season,  and  time.  The  six  climatological  regions  considered  are:  continental,  arctic, 
desert,  maritime,  tropical,  and  coastal.  The  three  time  periods  considered  are:  midnight,  sunrise,  and  local  noon. 

1 .3  Method.  In  this  study,  we  computed  probabilities  for  "generic"  climates  rather  than  for  specific  regions 
or  stations.  The  irregular  spacing  of  weather  stations,  differences  in  topography,  and  differences  in  instrument 
exposure  often  make  using  actual  data  difficult.  Additional  considerations  are  required  to  determine  if,  and 
to  what  extent,  the  observed  spatial  decorrelation  is  a  result  of  effects  common  to  similar  climates,  or  if  it  is 
a  result  of  local  effects.  To  make  this  determination,  we  used  a  statistical  technique  to  evaluate  the  spatial 
and  temporal  correlation  of  surface  observations. 

1.4  Output  Format.  The  probability  data  is  provided  on  a  series  of  charts  that  show  percentiles  of  the 
cumulative  distributions  of  changes  in  temperature,  wind  speed,  or  wind  direction  as  functions  of  distance  or 
time. 
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2.  STATISTICAL  MODEL 

2.1  Introduction.  The  evaluation  of  changes  in  temperature  or  wind  variables  as  a  function  of  distance  or 
time  is  complicated  by  the  irregular  spacing  of  weather  stations.  In  any  given  region,  station  spacing  is  roughly 
on  the  order  of  one  per  100  km,  which  provides  a  limited  sampling.  Although  there  are  databases  with  finer 
grid  spacing,  resolution  will  always  be  limited  by  station  spacing,  a  fact  that  makes  it  hard  to  analyze  the 
accuracy  of  the  data.  In  addition,  no  geographical  region  is  truly  homogeneous;  all  regions  have,  to  varying 
degrees,  local  features  that  affect  weather  observations.  Although  it’s  difficult  to  filter  out  these  microscale 
and  mesoscale  effects,  especially  for  winds,  we  devised  a  statistical  model  to  evaluate  changes  in  temperature, 
wind  speed,  and  wind  direction  with  distance.  We  selected  a  representative  station  for  each  climatological 
region.  Topeka,  Kansas,  for  example,  was  selected  to  represent  the  continental  climate.  In  this  model,  it  was 
assumed  that  all  points  surrounding  Topeka  have  exactly  the  same  climate;  for  each  model  climate,  the  mean 
and  standard  deviation  of  temperature  and  the  u-  and  v-components  of  the  wind  were  then  computed.  We 
assumed  that  these  variables  are  normally  distributed,  and  confirmed  this  assumption  with  tests  on  the 
skewness  and  kurtosis  of  the  actual  distributions. 

2.1 .1  Exception  for  Coastal  Climates.  We  cannot  make  the  assumption  of  homogeneity  in  coastal  areas 
because  of  discountinuities  resulting  from  the  sea-breeze  front,  which  is  largely  responsible  for  observed 
changes  in  temperature,  wind  speed,  and  wind  direction.  We  analyzed  the  coastal  climate  separately— see 
Chapter  4. 

2.2  Normal  Distributions.  The  assumption  that  temperatures  follow  a  normal  distribution  and  that  the  u- 
and  v-components  of  the  wind  follow  a  bivariate  normal  distribution  are  made  frequently  (Boehm,  1976;  Buell, 
1972;  Thiebaux  and  Pedder,  1987).  In  this  study,  we  assumed  that  the  joint  distribution  of  temperature 
between  two  points  (either  in  space  or  in  time)  is  bivariate  normal;  we  also  assumed  that  the  joint  distribution 
of  u-  and  v-components  of  wind  between  two  different  points  is  quadrivariate  normal.  The  fact  that  the 
marginal  distributions  are  normal  is  a  necessary  (but  not  sufficient)  condition  that  the  joint  probabilities  are 
also  normal  (Smith,  1976);  even  so,  this  assumption  is  often  made.  Univariate  normal  distributions  are  based 
on  two  parameters:  the  mean  and  standard  deviation.  The  bivariate  normal  distribution  has  two  means,  two 
standard  deviations,  and  a  fifth  parameter,  which  is  the  cross-correlation  coefficient.  Correlations  for  the 
quadrivariate  normal  distribution  are  more  complex  and  are  expressed  by  a  4x4  correlation  matrix  (see  Smith, 
1976,  for  further  discussion). 

2.3  Derivation.  Values  of  fluid  dynamical  variables  such  as  temperature  and  wind  speed  do  not  occur 
randomly;  instead,  they  exhibit  a  spatial  and  temporal  coherence  caused  by  stresses  within  the  viscous  fluid. 
One  would  expect,  for  example,  that  temperature  a  very  small  distance  away  from  temperature  7  at  a  central 
point  would  be  strongly  related  to  the  temperature  at  the  central  point.  In  other  words,  the  two  temperature 
measurements  should  be  strongly  positively  correlated  and  the  probability  of  a  large  temperature  difference 
|r  -  Tpl  between  the  two  points  should  be  small.  Correspondingly,  the  temperature  a  substantial  distance 
away  from  the  central  point  will  be  less  strongly  related  to  the  temperature  at  the  central  point;  the  two 
temperatures  are  generally  less  strongly  positively  correlated,  and  the  probability  of  a  large  temperature 
difference  between  the  two  is  correspondingly  larger.  The  key  point  is  that  the  correlation  p  between  T  and 
To  is,  in  general,  a  function  of  the  distance  r  between  the  two  points;  i.e.,  p  =  p(r).  The  temperature  is 
normally  distributed,  and  differences  between  temperatures  are  also  normally  distributed.  The  two  variables 
are  jointly  related  by  the  bivariate  normal  distribution,  with  one  random  variate  being  T,  the  other  being  T^, 
and  the  correlation  between  them  being  p(r).  If  the  temperature  at  a  given  point  is  T#,  the  probability  that 
a  temperature  7  at  a  distance  r  from  the  central  point  exceeds  Tq+AT  is  given  by  (Boehm,  1976): 
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(1) 


p(r  >  =  1  -  ^ 


-  Pir)p 

Vi  -  pW". 


where is  the  equivalent  normal  deviate  (END)  of  7,;  jl,  is  the  END  of  T,  +  AT;  and  ^(Q)  is  the  (univariate) 
normal  cumulative  probability  of  obtaining  an  END  of  less  than  Q  (see  Figure  1).  (For  a  discussion  of  the 
concept  of  equivalent  normal  deviates,  see  Whiton  and  Berecek,  1982.) 


Figure  1 .  Graphic  depiction  of  a  normai  probabiiity  density  function  and  the 
probabilities  described  by  Equations  1  and  2. 


It  must  be  emphasized  that  we  have  explicitly  stated  that  the  temperature  at  the  central  point,  T*,  has  been 
held  constant.  Various  values  of  T,  r,  and  AT  are  inserted  into  the  equation.  Next  we  compute  the  conditional 
probability  of  obtaining  a  temperature  less  than  T^-  AT: 


P{T<  Tf^-^T\T^  =  ^ 


y.  -  Pir)p 
Vi  -  P(rf. 


(2) 


where  is  the  END  of  To  -  AT.  Adding  these  cumulative  values  and  then  subtracting  this  sum  from  1,  yields 
the  conditional  probability  that  the  temperature  falls  within  a  range  of  2 AT  about  7#: 


P(\T  -Tol  <  2^T\T^  =  1  -  [P(T  >  To+^T\T^  +  P(T  <  Tq-^T\T^] 


(3) 
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To  find  the  unconditional  probability  that  the  difference  in  temperatures  between  the  two  points  is  less  than 
the  value  2 AT,  integrate  overall  values  of  temperature  at  the  central  point  (Tq): 


Pi\T  -  rj  <  2A7)  =  fP(\T  -rj  <  2AT\T^P(T^)dT^ 


(4) 


In  these  equations  we  account  for  the  effect  of  distance  through  the  correlation  coefficient  p=p(r).  It  is 
immaterial  how  one  subdivides  the  interval  of  temperature;  it  is  the  total  size  of  the  interval  that  matters.  We 
can  define  a  new  interval  AT'  =  2AT,  for  which  equations  3  and  4  will  be  valid. 

2.4  Wind  Speed  and  Direction.  Computations  of  the  changes  in  these  variables  with  distance  (or  time)  are 
complicated  by  the  required  variable  transformation.  Distributions  of  wind  speed  (V)  and  wind  direction  ( d) 
are  far  from  normal.  We  must  therefore  transform  V-0  components  to  «-v  components.  It  is  reasonable  to 
assume  that  the  joint  distribution  of  u  and  v  is  bivariate  normal  (Smith,  1976).  Computations  of  the  mean 
and  standard  deviations  of  u  and  v  and  the  cross  correlation  parameter  p  will  let  us  specify  the  probability 
density  function  (PDF),  p(u,v).  The  joint  probability  of  u  and  v  between  two  points  will  be  assumed  to  follow 
the  quadrivariate  normal  distribution.  There  are  two  basic  types  of  correlation  coefficients  in  this  type  of 
distribution.  The  first  is  the  autocorrelation  which  defines  the  correlation  between  the  two  u  (or  two  v) 
components  at  two  different  points  (either  in  space  or  time).  The  second  type  is  the  cross-correlation 
co^cient,  which  defines  the  correlation  between  the  u  and  v  components.  This  is  a  simplification,  since  the 
correlation  structure  in  the  quadrivariate  normal  distribution  is  a  4x4  matrix.  However,  following  the 
technique  in  Smith  (1976),  all  elements  of  this  matrix  can  be  obtained  by  specifying  the  two  auto-correlations: 

and  p,^  and  the  cross-correlation  p,^.  This  technique  ignores  the  so  called  "cross-lag  coefficients,"  which 
are  usually  small  and  difficult  to  estimate  from  real  data. 

2.4.1  Variable  Transformations.  Once  the  PDFs  for  the  components  of  the  wind  vector  have  been 
determined,  transformations  are  necessary  to  compute  the  corresponding  PDFs  for  wind  direction  and  speed. 
Conceptually,  this  can  be  carried  out  using  a  "checkerboard"  approach,  in  which  finite  increments  of  u  are 
shown  along  the  rows  and  finite  intervals  of  v  along  the  columns.  The  entry  in  each  checkerboard,  then,  is 
the  probability  of  the  joint  occurrence  of  the  given  u  and  v  intervals. 

2.4.2  Sample  Calculation.  To  show  how  this  technique  works,  consider  a  portion  of  the  checkerboard  as 

shown  in  Table  1.  Values  of  u  and  v  are  depicted  in  intervals  of  0.2  knots.  Suppose  the  probability  of 
(Mi=1.0,  Vi=0.8,  M2=1.2,  V2=1.0)  is  computed  to  be  0.02.  The  change  in  wind  direction  represented  by  these 
values  (from  Table  1)  is  11.5°.  We  use  this  information  to  compute  the  cumulative  distribution  function 
(CDF)  of  the  change  in  wind  direction.  The  value  P{u^  =  1.0,  Vj=  0.8,  =  1.2,  V2  =  1.0)  =  0.02  contributes 

a  value  0.02  to  the  CDF  for  the  value  P(Ad  <  11.5°).  A  summation  is  then  made  for  all  possible  entries,  and 
the  complete  CDF,  independent  of  the  initial  wind  direction  is  obtained.  The  probability  that  the  change  in 
wind  direction  40  is  less  than  a  threshold  value  400  is  given  mathematically  by: 

M  N  M  N  (5) 

p(S6  <  A9^  = 

i=i  j=i  *=1  /=! 

where  there  are  M  increments  of  u  and  N  increments  of  v.  The  change  in  wind  direction  between  («i,  v,)  and 
(«2,  V2)  is  represented  by  Ad'.  A  similar  computation  is  carried  through  for  wind  speeds.  This  highly  iterative 
process  is  accomplished  quickly  on  the  USAFETAC  mainframe  computer,  an  IBM  3090. 
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Table  1.  The  corresponding  change  in  wind  direction  (degrees)  for  the  stated 
vaiues  of  tv,,  u^,  v^,  (knots). 


u,  =  1.0  V,  =  0.8 
Ua 


Vz 

1.0 

1.2 

1.4 

1.6 

1.0 

6.3 

1.1 

3.1 

6.7 

1.2 

11.5 

6.3 

1.9 

1.8 

1.4 

15.8 

10.7 

6.3 

2.5 

1.6 

19.3 

14.5 

10.2 

6.3 

2.5  Temporal  Correlation.  A  technique  similar  to  that  described  in  Section  2.3  is  used  to  compute  the 
temporal  correlations.  The  same  distributions  are  assumed  to  apply.  For  the  spatial  correlations,  we  assumed 
that  the  mean  and  standard  deviation  of  the  variables  of  interest  were  the  same  at  both  points.  For  the 
temporal  case,  these  distributions  parameters  are  functions  of  time.  The  statistics  are  computed  based  on  ±3 
hours  of  the  mean.  The  actual  data  used  to  compute  the  statistical  parameters  are  all  obtained  from  stations 
that  report  hourly.  Thus,  the  mean  and  standard  deviation  for  each  of  the  hours  of  interest  are  computed. 

2.6  Correlation  Coefficients.  The  spatial  correlation  of  temperature  as  a  function  of  distance  was  obtained 
from  Bertoni  and  Lund  (1964);  it  represents  the  average  of  many  hundreds  of  pairs  of  stations  in  Europe. 
Bertoni  and  Lund  obtained  a  funaional  relationship  of  correlation  (pj  with  distance  (r): 

p  (r)  =  exp(-— )  (<>) 

r. 


where  r,  is  the  scale  distance,  equal  to  1,000  km.  The  spatial  correlation  of  «-  and  v-components  were 
obtained  from  Buell  (1972),  who  suggested  the  model: 


p/r)  =  exp(-r2/r^) 


(7) 


where  the  scale  distance  r,  is  250  km.  Cross-correlations  and  temporal  correlations  were  all  computed  from 
observations  from  each  of  the  representative  stations-see  Section  2.8.  The  differences  in  the  corresponding 
scale  distance  reveal  that  the  wind  component  spatial  correlations  decrease  faster  with  distance  than  with 
temperature. 

2.6.1  Sensitivity  to  Correlation  Coefficient.  Clearly,  our  assumptions  of  a  universal  spatial  airrelation 
structure  are  significant.  Studies  of  the  changes  in  this  structure  with  various  geographical  regions,  or  as  a 
function  of  season,  have  not  been  made;  however,  high  precision  of  the  correlation  coefficient  is  not  required, 
as  will  be  shown  by  the  following  analysis. 

2.6.2  Meaning  of  the  Correlation  Coefficient.  For  any  normal  distribution,  68%  of  all  observations  lie 
within  one  standard  deviation  ( a)  of  the  mean.  Similarly,  50%  of  all  observations  lie  within  0.675  er  of  the 
mean.  For  two  distributions  with  the  same  mean  and  the  same  standard  deviation,  the  correlation  coefficient 
is  a  measure  of  the  reduction  in  the  variance  in  one  distribution  when  the  value  of  the  other  variable  is  known. 


5 


For  example,  consider  a  temperature  distribution  with  a  mean  of  90°  F  and  a  standard  deviation  of  6  F.  We 
will  assume  that  this  distribution  describes  the  temperature  at  two  points,  A  and  B.  If  the  temperature  at 
point  A  is  90°  F  and  the  autocorrelation  (p)  is  0.8,  the  standard  deviation  of  the  conditional  distribution  (<7,) 
of  temperature  at  point  B  is  given  by  (Smith,  1976): 

a.  =  o/T? 


or  in  this  example,  3.6°  F.  An  approximation  of  the  range  where  50%  of  the  temperatures  at  point  B  lie  is 
then  2.2°  F.  If  the  autocorrelation  were  as  low  as  0.6  the  range  would  be  3.2°  F.  Thus,  a  significant  change 
in  correlation  (0.8  versus  0.6)  results  in  only  a  1°  F  change  in  the  estimate  of  the  50th  percentile  o 
temperature  differences. 

2.7  Statistical  Parameters.  Although  this  example  shows  this  problem’s  limited  sensitivity  to  the  correlation 
coefficient,  it  also  highlights  the  much  larger  sensitivity  of  the  final  results  to  the  standard  deviation.  A  20% 
increase  in  standard  deviation  results  in  a  corresponding  20%  increase  in  the  range  of  temperatures  expected. 
Less  obvious  in  this  acample,  is  the  sensitivity  of  the  means  of  the  two  distributions.  For  example,  when 
evaluating  the  change  in  temperature  over  a  period  of  2  hours,  the  mean  temperature  at  each  hour  may 
sicnificantly  different.  Thus,  even  when  correlations  are  high  and  standard  deviations  are  low,  if  the  change 
in  the  mean  temperature  is  3°  F,  the  most  likely  difference  will  be  at  least  3°  F.  For  most  ^e  spatia 
calculations,  this  effect  was  not  important  because  we  assumed  the  means  were  constant.  For  temporal 
calculations,  however,  the  effect  was  significant. 

2  8  Station  Selection.  Six  stations  were  selected  for  this  study,  one  for  each  of  the  climatologiral  regions 
introduced  in  Section  1.2.  All  were  in  the  U.S.,  primarily  because  of  the  availability  of  ftoiirfy  wither 
observations  (at  many  overseas  locations,  weather  observations  are  available  only  every  3  hours-thts  makes 
determination  of  temporal  correlation  nearly  impossible).  The  U.S.  stations  were 

reliable-errors  in  the  measurement  of  temperature,  wind  speed,  and  wind  direction  (especially  biased  errors) 
are  extremely  difficult  to  detect.  A  list  of  the  stations  selected  is  shown  in  Table  2. 

Table  2.  Stations  used  as  representative  of  the  climates  selected  in  this  study. _ 


Block 

Climate 

Station 

Station  # 

Continental 

Forbes  Field  KS 

724565 

Arctic 

Eielson  AFB  AK 

702650 

Tropical 

Howard  AFB  PN 

788060 

Maritime 

Hickam  AFB  HI 

911803 

Desert 

Nellis  AFB  NV 

723865 

Coastal 

Patrick  AFB  FL 

747950 

Latitude 

38°  57’  N 
64°  40’  N 
8°  55’  N 
21°  19’  N 
36°  14’  N 
28°  14’  N 


Longitude  Elev(m) 


95°  40’  W 
147°  06’  W 
79°  36’  W 
157°44’ W 
115°  02’ W 
80°  36’  W 


329 

167 

16 

4 

570 

3 


Observations  from  these  stations  were  used  to  compute  the  mean  and  standard  deviations  of  temperature,  the 
u  and  V  wind  components,  the  wind  component  cross-correlations,  and  the  serial  correlations  for  temperature 

and  wind  components. 
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3.  VERIFICATION  OF  ASSUMPTIONS 

3.1  Test  for  Normality.  A  fundamental  assumption  of  the  model  is  that  temperature  and  the  u-  and  v- 
components  of  wind  are  all  normally  distributed.  To  test  for  this  assumption,  wc  used  a  technique  by 
D’Agostino  et  al.  (1990),  which  uses  observed  skewness  and  kurtosis  to  determine  whether  or  not  the  normal 
distribution  can  be  considered  an  approximation  to  the  observed  distribution.  The  chi-square  and  similar  tests 
are  much  stricter;  they  consider  the  observed  distribution  non-normal  if  any  type  of  statistically  significant 
difference  is  observed. 


3.1 .1  Skewness.  Skewness  is  a  measure  of  the  symmetry  of  the  distribution.  It  is  represented  mathematically 
as  bi  and  defined  as: 


(9) 


where 


m~ 


I 

n 


-  5)3 


(10) 


and 

=  (11) 

n 

In  these  equations,  Xj  represents  the  individual  observations,  x  is  the  mean  value,  and  n  is  the  number  of 
observations.  A  skewness  of  zero  represents  a  symmetrical  distribution. 


3.1.2  Kurtosis.  Although  kurtosis  is  difficult  to  define,  it  basically  discriminates  between  "peaked"  and  "flat- 
topped"  distributions.  It  is  assigned  the  symbol  bj  defined  as: 


where 


(12) 


m 


n  ,=i 


(13) 


3.1.3  interpretation  of  Skewness  and  Kurtosis.  By  definition,  a  normal  distribution  has  neither  skewness 
nor  kurtosis;  i.e.,  =  0  and  b2  =  0.  When  the  skewness  is  non-zero,  the  mean  and  the  median  do  not 

coincide.  When  the  kurtosis  is  positive,  the  distribution  is  flat  and  the  number  of  observations  that  occur 
within  one  standard  deviation  of  the  mean  is  too  small  for  a  normal  distribution.  Conversely,  a  negative 
kurtosis  reveals  that  too  many  observations  occur  within  one  standard  deviation  of  the  mean.  For  a  normal 
distribution  to  be  considered  a  reasonable  approximation  to  an  observed  distribution,  both  the  skewness  and 
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kurtosis  values  must  be  reasonably  close  to  zero.  The  criteria  presented  by  D’Agostino  et  al.  (1990)  to  test 
for  normality  are; 


IM  <  4 


and 


\h\  <  4 


When  both  conditions  are  met,  the  observed  distribution  can  be  reasonably  approximated  as  normal. 

3.2  Tests  on  Temperature.  Skewness  and  kurtosis  values  were  computed  for  observations  of  temperature 
at  Topeka  for  the  four  seasons  and  three  time  periods  used  in  this  study;  the  results  are  shown  in  Table  3. 
Also  listed  are  the  skewness  and  kurtosis  thresholds  (THSH)  for  testing  normality.  Skewness  was  detected 
at  two  time  periods:  spring-noon  and  summer-noon,  but  kurtosis  was  detected  only  at  summer-noon.  For 
these  time  periods,  the  assumption  of  normality  is  not  supported.  The  data  for  summer  noon  revealed  much 
larger  values  of  skewness  and  kurtosis  than  for  the  other  time  groups.  However,  it  generally  justifies  the 
assumption  that  temperature  is  normaliy  distributed. 


^  (15) 

n 


6  (14) 

n 


Table  3.  Skewness  and  kurtosis  values  for  temperature  distributions  at  Topeka.  Also  listed  are  the 
number  of  observations  (N)  and  thresholds  for  assuming  a  normal  distribution. _ 


Skew 

Kurt 

Criteria 

Season 

Time 

Skew 

Kurt 

N 

Thsh 

Thsh 

Met? 

WINTER 

Midnight 

-0.22 

-0.22 

725 

0.36 

0.73 

Y 

Sunrise 

-0.24 

-0.51 

835 

0.34 

0.68 

Y 

Noon 

0.06 

-0.42 

833 

0.34 

0.68 

Y 

SPRING 

Midnight 

0.13 

-0.42 

690 

0.37 

0.75 

Y 

Sunrise 

-0.17 

-0.46 

809 

0.34 

0.69 

Y 

Noon 

-0.53 

0.53 

810 

0.34 

0.69 

N 

SUMMER 

Midnight 

-0.26 

0.03 

688 

0.37 

0.75 

Y 

Sunrise 

-0.21 

-0.37 

811 

0.34 

0.69 

Y 

Noon 

-0.92 

1.76 

811 

0.34 

0.69 

N 

AUTUMN 

Midnight 

-0.02 

-0.59 

713 

0.37 

0.73 

Y 

Sunrise 

-0.03 

-0.42 

835 

0.34 

0.68 

Y 

Noon 

0.01 

-0.39 

836 

0.34 

0.68 

Y 

3.3  Tests  on  Wind  Components.  Similar  tests  were  carried  out  to  evaluate  the  assumption  of  normality  for 
the  u-  and  v-components  of  wind.  These  results  are  shown  in  Tables  4  and  5,  respectively.  It’s  clear  that  there 
are  times  when  the  assumption  of  normality  is  not  justified.  In  fact,  the  test  fails  for  the  u-component  during 
all  four  seasons  at  midnight.  The  large  positive  values  of  kurtosis  indicate  a  peaked  distribution;  too  many 
observations  are  clustered  around  the  mean.  This  results  from  the  fact  that  calm  winds  occur  with  higher 
frequency  at  night.  Calm  (and  light)  winds  create  additional  problems.  The  transformation  between  u-v 
coordinates  and  V-6  coordinates  can  become  unstable  with  light  winds.  A  small  change  in  m  or  v  can  result 
in  large  changes  in  6.  EXAMPLE:  The  components  («  =  1.0,  v  =  0.5)  corresponds  to  (0  =  243°,  V  =  1.1). 
A  slight  increase  in  the  v-component  to  1.5  knots  results  in  a  wind  direction  of  213°  and  a  wind  speed  of  1.8 
knots;  a  1-knot  change  in  one  of  the  components  resulted  in  a  30°  change  in  wind  direction. 

To  solve  these  problems,  the  statistics  were  recomputed  to  exclude  wind  speeds  below  3  knots.  This  reduced 
kurtosis  values,  but  generally  not  below  the  threshold  values.  It  also  resulted  in  an  increase  in  kurtosis  values 
for  the  v-components.  For  example,  the  kurtosis  of  the  w-component  at  winter-midnight  decreased  from  2.14 
to  1.36.  However,  the  kurtosis  of  the  v-component  changed  from  0.12  to  -0.42.  The  assumption  of  normality 
is  clearly  a  problem  for  modeling  light  winds.  Since  light  winds  are  not  of  much  interest  in  this  study, 
computed  statistics  exclude  conditions  when  winds  were  less  than  3  knots.  It  is  clear  that  the  results  of  this 
study  do  not  apply  when  winds  are  light. 


Table  4.  Values  of  skewness  and  kurtosis  for  the  u-component  of  surface  wind  at  Topeka.  Also  listed 
are  the  number  of  obsen/ations  (N)  and  thresholds  for  assuming  a  normal  distribution. 


Season 

Time 

Skew 

WINTER 

Midnight 

0.26 

Spring 

-0.01 

Noon 

0.27 

SPRING 

Midnight 

0.18 

Spring 

0.13 

Noon 

0.16 

SUMMER 

Midnight 

-0.26 

Spring 

-0.15 

Noon 

0.35 

AUTUMN 

Midnight 

0.31 

Spring 

0.19 

Noon 

0.52 

Kurt 

JL 

Skew 

Thsh 

Kurt 

Thsh 

Criteria 

Met? 

2.14 

723 

0.36 

0.73 

N 

0.36 

836 

0.34 

0.68 

Y 

0.46 

836 

0.34 

0.68 

Y 

1.45 

690 

0.37 

0.75 

Y 

0.15 

810 

0.34 

0.69 

Y 

0.20 

810 

0.34 

0.69 

Y 

2.83 

713 

0.37 

0.73 

N 

0.41 

836 

0.34 

0.68 

Y 

1.58 

837 

0.34 

0.68 

N 

2.65 

713 

0.37 

0.73 

N 

0.65 

837 

0.34 

0.68 

Y 

1.14 

837 

0.34 

0.68 

N 
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Table  5.  Values  of  the  skewness  and  kurtosis  of  the  v'-component  of  wind  at  Topeka.  Also  listed  are 
the  number  of  observations  (N)  and  thesholds  for  assuming  a  normal  distribution. 


Skew 

Kurt 

Criteria 

Season 

Time 

Skew 

Kurt 

N 

Thsh 

Thsh 

Met? 

WINTER 

Midnight 

0.00 

0.12 

723 

0.36 

0.73 

Y 

Sunrise 

-0.01 

-0.40 

836 

0.34 

0.68 

Y 

Noon 

0.10 

-0.55 

836 

0.34 

0.68 

Y 

SPRING 

Midnight 

0.12 

0.47 

690 

0.37 

0.75 

Y 

Sunrise 

0.19 

-0.81 

810 

0.34 

0.69 

N 

Noon 

0.30 

-0.68 

810 

0.34 

0.69 

Y 

SUMMER 

Midnight 

0.27 

0.82 

713 

0.37 

0.73 

N 

Sunrise 

-0.14 

-0.27 

836 

0.34 

0.68 

Y 

Noon 

-0.14 

-0.50 

837 

0.34 

0.68 

Y 

AUTUMN 

Midnight 

0.23 

0.73 

713 

0.37 

0.73 

Y 

Sunrise 

-0.05 

-0.55 

837 

0.34 

0.68 

Y 

Noon 

0.01 

-0.66 

837 

0.34 

0.68 

Y 

3.4  MuKivariate  Normal  Distributions.  To  validate  the  assumptions  that  the  temperature  and  u-  and  v- 
componcnts  of  wind  at  two  points  arc  distributed  according  to  a  multivariate  normal  distribution,  sample 
model  distributions  and  an  empirically  computed  distribution  were  compared.  Figure  2  shows  a  comparison 
of  the  cumulative  distributions  of  the  3-hour  change  in  temperature  during  the  winter-sunrise  period  at  Topeka 
with  the  3-hour  change  from  the  statistical  model.  The  winter-sunrise  period  was  selected  because  it  had  the 
largest  amount  of  variability  of  all  12  time  periods  studied.  Figure  3  shows  the  3-hour  change  in  wind 
direction  at  the  same  time  period,  and  Figure  4  shows  the  3-hour  change  in  wind  speed.  Similar  graphics 
(Figures  5-7)  were  prepared  for  the  spatial  correlation. 


0  t  4eaioi2i4 
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Figure  2.  Comparison  of  the  observed  and 
simulated  cumulative  distribution  functions  for 
the  3-hour  change  in  temperature  for  Topeka 
at  winter-sunrise. 
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Figure  3.  Comparison  of  the  observed  and 
simulated  cumulative  distribution  functions  for 
the  3-hour  change  in  wind  direction  for 
Topeka  at  winter-sunrise. 
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Figure  4.  Comparison  of  the  observed  and 
simulated  cumulative  distribution  functions  for 
the  3-hour  change  in  wind  speed  for  Topeka, 
at  winter-sunrise. 


Figure  5.  Comparison  of  the  observed  and 
simulated  cumulative  distribution  functions  of 
the  difference  in  temperature  between  Topeka 
and  Ft  Riley  at  winter-sunrise. 
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Figure  6.  Comparison  of  the  observed  and 
simulated  cumulative  distribution  functions  for 
the  difference  in  wind  direction  between 
Topeka  and  Ft  Riley  at  winter-sunrise. 
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Figure  7.  Comparison  of  the  observed  and 
simulated  cumulative  distribution  functions  for 
the  difference  in  wind  speed  between  Topeka 
and  Ft  Riley  at  winter-sunrise. 


The  changes  in  temperature,  wind  direction,  and  wind  speed  between  Topeka  and  Fort  Riley  (60  nm  west  of 
Topeka)  were  computed  and  compared  against  the  statistical  model;  there  was  good  agreement.  Model 
estimates  of  the  percentile  values  through  about  the  75th  percentile  are  very  close  to  the  observed  values. 
However,  estimates  above  the  75th  percentile  begin  to  deviate,  due  mainly  to  the  fact  that  we  are  now  using 
values  in  the  tails  of  the  modeled  distributions.  Departures  from  normality  are  most  notable  in  these  tails, 
and  small  errors  are  amplified.  This  effect  is  further  enhanced  by  the  small  slope  in  the  curve  beyond  about 
the  75th  percentile.  EXAMPLE:  The  model  estimate  of  the  99th  percentile  of  the  3-hour  change  in 
temperature  is  14.3  and  the  actual  value  is  11.9.  This  may  seem  to  be  a  large  error,  but  the  model  estimate 
for  the  corresponding  percentile  to  11.9  is  98.5.  In  this  region,  small  changes  in  percentile  result  in  large 
changes  in  temperature.  We  do  not  believe  that  this  effect  is  significant.  Estimates  of  values  even  out  to  the 
95th  percentile  are  reasonable. 
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4.  COASTAL  ENVIRONMENT 


4.1  Climate.  The  assumption  of  large-scale  climatic  homogeneity  is  clearly  not  true  in  coastal  environments 
because  the  effect  of  the  sea-breeze  front  provides  a  sharp  discontinuity  in  air  mass  over  relatively  short 
distances.  The  movement  of  this  front  adds  additional  variability  to  atastal  climates.  We  will  attempt  to 
consider  the  change  in  climate  with  distance  to  evaluate  changes  in  temperature,  wind  speed,  and  wind 
direction  with  distance. 


4.2  Mesoscale  Correlation  Structure.  Since 
large  temperature  changes  occur  over  relatively 
short  distances,  we  need  to  be  familiar  with  the 
correlation  structure  in  the  mesoscale.  To  obtain 
this  information,  we  used  the  Weather  Information 
Network  Display  System  (WINDS)  dataset  from 
Cape  Canaveral,  Florida.  The  WINDS  network 
consists  of  28  sensors  that  measure  temperature, 
wind  direction,  and  speed.  The  largest  distance 
between  any  two  sensors  is  30  km.  Figure  8  is  a 
map  of  sensor  locations.  Temperature  and  wind 
data  are  available  for  various  levels,  but  for  this 
study  we  considered  only  the  lowest  (10-meter) 
level.  Next,  we  tried  to  find  a  relationship  of 
temperature  correlation  with  distance.  A  matrix  of 
temperature  correlations  was  computed  for  all 
possible  pairs  of  stations.  This  can  be  depicted  as: 


Figure  8.  Map  of  Patrick  AFB  WINDS  sensors. 


P  = 


1  Pl2  Pl3  Pl4  Pl5 

1  P23  P24  P25 

1  P34  P35 

•  •  •  1  P45 

.  .  .  .  1 


(16) 


where  p  is  the  matrix  of  correlation  coefficients.  The  members  of  the  matrix  represent  the  correlation 
between  two  locations.  For  example,  the  correlation  between  locations  3  and  5  is  given  by  P35.  Only  the  terms 
above  the  diagonal  need  to  be  computed,  since  this  matrix  is  symmetric  about  the  diagonal  (p,^  =  p^,,  etc.). 
Corresponding  to  each  correlation,  (e.g.,  P35),  is  a  distance,  (e.g.,  rjj),  between  the  two  sensors.  This  was 
carried  out  for  all  the  WINDS  stations,  llie  correlations  are  based  on  hourly  temperature  observations. 
Figure  9  is  a  plot  of  these  correlations  of  temperature  obtained  from  the  WINDS  dataset  as  a  function  of 
distance.  Also  shown  is  the  expected  correlation  structure  using  the  model  of  Bertoni  and  Lund  (1964)~see 
equation  6.  A  functional  relationship  between  correlation  and  distance  is  not  observed  in  the  mesoscale. 
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Figure  9.  Correlation  of  temperature  among 
Patrick  AFB  WINDS  locations  as  a  function  of 
distance.  Also  shown  is  the  expected  correlation 
using  the  model  of  Bertoni  and  Lund. 


Figure  10.  Correlation  of  temperature  between 
Patrick  AFB  and  neighboring  synoptic  stations. 


A  regression  fit  of  these  points  resulted  in  a  fine  independent  of  distance,  suggesting  that  distance  alone  does 
not  describe  the  decorrelation  of  temperature  observations  at  the  mesoscale.  Differences  caused  by  orographic 
controls  and  the  sea  breeze  appear  to  be  much  more  significant  than  distance  once  the  separation  becomes 
less  than  critical.  We  propose  that  this  critical  distance  be  used  to  define  the  mesoscale.  Similar  behavior 
is  observed  in  the  wind  data.  Other  information  is  required  to  compute  the  correlation  coefficient  within  the 
mesoscale.  Tabor  and  Motte  (1988)  found  that  a  detailed  terrain  dataset  was  sufficient  to  extrapolate  wind 
speed  and  direction.  This  suggests  that  a  mesoscale  analysis  of  weather  variables  independent  of  a  specific 
location  is  not  practical. 

4.3  Larger-Scale  Correlation  Structure.  To  supplement  the  mesoscale  network,  we  used  synoptic  scale 
stations  to  determine  whether  the  correlation  models  described  in  equations  6  and  7  also  applied  in  coastal 
environments.  The  results  of  this  study,  which  showed  that  the  models  were  indeed  valid,  are  shown  in  Figure 
10.  The  scale  distance  in  this  plot  (about  1,300  km)  is  slightly  larger  than  the  value  (1,000  km)  given  in 
Section  2.7.  This  larger  scale  distance  was  used  in  carrying  out  the  calculations  for  the  coastal  environment. 
The  wind  structure  showed  similar  behavior,  but  the  data  from  the  limited  number  of  stations  was  not  as 
smooth.  Therefore,  we  used  the  same  scale  distance  as  described  in  Section  2.7  (250  km)  for  the  wind 
calculations. 


4.4  Changes  in  Mean  Temperature.  As  shown  in  Figures  11-14,  the  north-south  and  east-west  temperature 
gradients  are  both  significant.  Figure  11  shows  January  mean  maximum  temperatures  throughout  Florida; 
Figure  12,  mean  minimums.  Figure  13  shows  July  mean  maximum  temperatures;  Figure  14,  mean  minimums. 
These  figures  show  the  effects  of  the  ocean,  which  are  much  more  pronounced  in  winter  than  in  summer. 
Maximum  temperatures  during  the  day  are  warmer  inland,  but  at  night  the  warmest  temperatures  are  along 
the  coast.  There  was  little  change  in  the  standard  deviation  of  the  temperature  distributions  at  any  given  hour 
as  a  function  of  distance.  These  figures  are  based  on  data  in  Bradley  (1978). 

4.5  Distribution  of  Mean  Temperature.  Figures  11-14  show  the  observed  change  in  mean  temperature  as 
a  function  of  distance  inland.  Again,  the  assumption  of  a  uniform  mean  (made  for  the  other  climates)  clearly 
does  not  apply.  We  therefore  assume  the  mean  to  be  a  function  of  distance-,  an  assumption  of  a  constant 
standard  deviation  is  reasonable.  We  will  modify  our  statistical  model  by  including  a  relationship  for  the 
change  in  mean  temperature  as  a  function  of  distance  in  coastal  climates. 
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Figure  11.  Mean  maximum  temperature  Figure  12.  Mean  minimum  temperature 
throughout  Fiorida  in  January.  throughout  Florida  in  January. 


Figure  13.  Mean  maximum  temperature  temperature 

throughout  Fiorida  in  July.  throughout  Florida  in  July. 


Figure  15  shows  the  modeled  change  in  mean  summer  maximum  temperature  as  one  moves  inland.  We 
assumed  a  constant  slope  out  to  a  distance  of  80  km,  and  a  constant  value  beyond  that.  The  data  in  Figure 
15,  based  on  the  maximum  temperature  plots,  was  applied  to  calculations  for  changes  in  temperature  at  noon. 
Figure  16  is  similar,  but  based  on  minimum  temperatures;  this  model  was  applied  to  calculations  made  for 
sunrise.  We  assumed  a  constant  mean  for  midnight,  estimating  that  this  was  about  the  midpoint  between 
maximum  and  minimum  temperatures.  Figures  17  and  18  depict  the  corrections  applied  for  winter.  For  spring 
and  autumn,  averages  of  the  corrections  presented  for  winter  and  summer  were  used. 


Figure  15.  Modeled  change  in  mean  Figure  16.  Modeled  change  in  mean 

temperature  as  a  function  of  distance  inland  temperature  as  a  function  of  distance  inland 

during  summer  at  noon.  during  summer  at  sunrise. 


Figure  17.  Modeled  change  in  mean 
temperature  as  a  function  of  distance  inland 
during  winter  at  noon. 


Figure  18.  Modeled  change  in  mean 
temperature  as  a  function  of  distance  inland 
during  winter  at  sunrise. 


4.6  Distribution  of  the  u-  and  v-Components  of  Wind.  Our  analysis  of  mean  u-  and  v-components  as  a 
function  of  distance  did  not  show  a  pattern  similar  to  temperature.  Since  there  was  no  functional  change  in 
the  mean  wind  components  as  a  function  of  distance,  we  have  assumed  constant  values  for  these  variables. 

4.7  Conclusion.  The  statistical  model  that  describes  the  relationship  of  temperature  (or  wind  direction  or 
wind  speed)  cannot  be  used  to  describe  mesoscale  relationships.  At  scales  less  than  some  critical  distance, 
correlation  becomes  independent  of  distance  and  dependent  on  the  local  orographic  controls  as  they  interact 
with  the  atmosphere  to  produce  mesoscale  effects.  To  describe  mesoscale  relationships,  a  mesoscale  model 
must  be  used.  The  results  of  such  a  model  is  restricted  to  the  particular  location  being  studied.  A  generic 
approach  to  mesoscale  analysis  (not  tied  to  any  particular  location)  is  not  practical. 
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5.  MODEL  OUTPUT 


5.1  Results.  Model  output  is  given  as  figures  in  Appendices  A-F.  For  each  climate,  there  are  18  figures 
per  >veather  element.  Each  page  depicts  one  of  three  time  periods  (local  midnight,  sunrise,  and  local  noon), 
for  each  of  the  four  seasons.  Seasonal  statistics  were  based  on  a  representative  month  (winter  =  January; 
spring  =  April;  summer  =  July;  autumn  =  October).  Appendices  A-C  give  spatial  results  (A-temperature, 
B-wind  direction,  C-wind  speed).  Appendices  D-F  give  temporal  results  (D-temperature,  E-wind  direction,  F- 
wind  speed).  Altogether,  there  are  432  graphs  (four  per  page;  18  pages  per  appendix;  six  appendices). 
Although  the  charts  are  labeled  as  correlations,  there  are  no  true  correlation  statistics  in  the  graphs. 

5.2  Interpretation.  The  horizontal  axis  in  the  graphs  depicting  temporal  results  gives  the  change  in  time 
from  -3  hours  to  4-3  hours.  The  vertical  axis  gives  the  change  in  each  variable  of  interest  (temperature,  wind 
direction,  or  wind  speed).  For  the  spatial  correlations,  the  horizontal  axis  lists  the  change  in  distance  from 
0  to  250  km;  the  vertical  axis  is  the  same.  There  are  four  graphs  on  each  page.  The  climatological  region  and 
time  period  are  constant,  but  each  graph  depicts  a  different  season.  On  each  graph,  five  curves  are  plotted: 
from  bottom  to  top,  they  are  the  25th,  50th,  75th,  90th,  and  95th  percentiles  of  y-axis  values. 

5.2.1  Example  In  the  Temporal  Case.  Refer  to  the  first  page  in  Appendix  D,  and  to  the  upper  left-hand 
plot  (the  temporal  correlation  of  temperature  for  continental  climates  at  winter/midnight).  The  estimates  for 
the  percentiles  of  the  synoptic-scale  change  in  temperature  from  3  hours  before  midnight  to  midnight  are: 
25th  percentile~0.5°  C;  50th  percentile- 1.0°  C;  75th  percentile-2.0°  C;  90th  percentile-2.9°  C;  and  95th 
percentile-4.5°  C.  In  other  words,  at  continental  winter  midnight,  considering  synoptic-scale  processes  alone, 
there  is  only  a  10%  risk  of  encountering  a  3-hour  temperature  change  greater  than  2.9°  C  and  only  a  5%  risk 
of  a  change  greater  than  4.5°  C.  If  that  risk  is  unacceptably  high,  an  observing  interval  more  frequent  than 
once  every  3  hours  is  recommended. 

5.2.2  Example  In  the  Spatial  Case.  Consider  Appendbc  A,  which  shows  continental  winter  midnight 
temperature  differences  as  a  function  of  the  distance  between  two  observing  stations.  Within  50  km, 
considering  synoptic-scale  processes  alone,  there  is  about  an  18%  risk  of  encountering  a  temperature  difference 
of  3°  C,  whereas  at  75  km  there  is  a  25%  risk  of  a  temperature  difference  that  large,  and  at  100  km  there  is 
a  50-50  chance  of  encountering  such  a  difference.  If  a  system  under  design  requires  that  a  field  of 
temperatures  be  observed  and  that  90%  of  the  observations  must  be  within  an  error  of  3°  C  (not  considering 
mesoscale  features),  then  the  observing  network  must  be  spaced  no  wider  than  about  37  km.  At  continental 
summer  midnight,  the  same  performance  can  be  achieved  with  a  network  spacing  of  175  km. 

5.2.3  Additional  Example.  If  at  continental  winter  midnight  it  is  necessary  to  estimate  the  temperature  at 
a  point  where  there  is  no  weather  station  by  using  data  from  a  point  200  km  away,  then  the  median  error, 
considering  synoptic-scale  processes  alone  is  3°  C,  and  90%  of  such  errors  will  be  within  7°  C.  At  continental 
summer  midnight,  the  median  error  in  using  a  station  200-km  distant  is  well  within  instrument  observing  error, 
and  90%  of  such  errors  are  within  3°  C. 

5.2.4  Plots  with  Little  Variability.  In  some  of  the  charts,  the  lines  are  very  close  together.  This  is 
particularly  true  for  plots  of  the  maritime  climate,  in  which  the  differences  between  any  of  the  percentiles  are 
small.  The  top  curve  always  represents  estimates  of  the  95th  percentile. 

5.3  Final  Point.  Although  this  study  was  limited  to  the  spatial  and  temporal  correlation  of  only  three  surface 
weather  variables,  the  technique  can  be  applied  to  other  variables,  both  surface  and  upper-air.  It  is  especially 
suited  for  upper-air  variables  where  mesoscale  effects  are  less  significant.  Smith  (1976)  used  this  technique 
to  interpolate  wind  observations  vertically,  suggesting  that  a  study  on  the  vertical  correlation  of  other  weather 
variables  might  also  be  possible. 
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GLOSSARY 


CDF 

cumulative  distribution  function 

END 

equivalent  normal  deviate 

km 

kilometer 

NM 

nautical  mile 

PD 

probability  distribution  function 

Thsh 

Threshold 

USAFETAC 

United  States  Air  Force  Environmental  Technical  Applications  Center 

WINDS 

Weather  Information  Network  Display  System 

19 


APPENDIX  A 


SPATIAL  CORRELATION  OF  TEMPERATURE 

A-1  Spatial  correlation  of  temperature,  continental  climate,  midnight. 
A-2  Spatial  correlation  of  temperature,  continental  climate,  sunrise. 
A-3  Spatial  correlation  of  temperature,  continental  climate,  noon. 

A-4  Spatial  correlation  of  temperature,  arctic  climate,  midnight. 

A-5  Spatial  correlation  of  temperature,  arctic  climate,  sunrise. 

A-6  Spatial  correlation  of  temperature,  arctic  climate,  noon. 

A-7  Spatial  correlation  of  temperature,  desert  climate,  midnight. 

A-8  Spatial  correlation  of  temperature,  desert  climate,  sunrise. 

A-9  Spatial  correlation  of  temperature,  desert  climate,  noon. 

A-10  Spatial  correlation  of  temperature,  maritime  climate,  midnight. 
A-1 1  Spatial  correlation  of  temperature,  maritime  climate,  sunrise. 
A-1 2  Spatial  correlation  of  temperature,  maritime  climate,  noon. 

A-1 3  Spatial  correlation  of  temperature,  tropical  climate,  midnight. 
A-1 4  Spatial  correlation  of  temperature,  tropical  climate,  sunrise. 

A-1 5  Spatial  correlation  of  temperature,  tropical  climate,  noon. 

A-1 6  Spatial  correlation  of  temperature,  coastal  climate,  midnight. 
A-1 7  Spatial  correlation  of  temperature,  coastal  climate,  sunrise. 

A-1 8  Spatial  correlation  of  temperature,  coastal  climate,  noon. 
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A-8  Spatial  correlation  of  temperature,  desert  climate,  sunrise. 
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A-9  Spatial  correlation  of  temperature,  desert  climate,  noon. 
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A-10  Spatial  correlation  of  temperature,  maritime  climate,  midnight. 
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A-1 1  Spatial  correlation  of  temperature,  maritime  climate,  sunrise. 
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A-14  Spatial  correlation  of  temperature,  tropical  climate,  sunrise. 
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A-15  Spatial  correlation  of  temperature,  tropical  climate,  noon. 


A-16  Spatial  correlation  of  temperature,  coastal  climate,  midnight. 
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A-17  Spatial  correlation  of  temperature,  coastal  climate,  sunrise. 


APPENDIX  B 


SPATIAL  CORRELATION  OF  WIND  DIRECTION 

B-1  Spatial  correlation  of  wind  direction,  continentai  climate,  midnight. 
B-2  Spatial  correlation  of  wind  direction,  continentai  climate,  sunrise. 
B-3  Spatial  correlation  of  wind  direction,  continental  climate,  noon. 

B-4  Spatial  correlation  of  wind  direction,  arctic  climate,  midnight. 

B-5  Spatial  correlation  of  wind  direction,  arctic  climate,  sunrise. 

B-6  Spatial  correlation  of  wind  direction,  arctic  climate,  noon. 

B-7  Spatial  correlation  of  wind  direction,  desert  ciimate,  midnight. 

B-8  Spatial  correlation  of  wind  direction,  desert  climate,  sunrise. 

B-9  Spatial  correlation  of  wind  direction,  desert  ciimate,  noon. 

B-10  Spatial  correlation  of  wind  direction,  maritime  climate,  midnight. 
B-1 1  Spatial  correlation  of  wind  direction,  maritime  climate,  sunrise. 
B-1 2  Spatial  correlation  of  wind  direction,  maritime  climate,  noon. 

B-1 3  Spatial  correlation  of  wind  direction,  tropicai  climate,  midnight. 
B-1 4  Spatial  correlation  of  wind  direction,  tropical  climate,  sunrise. 

B-1 5  Spatial  correlation  of  wind  direction,  tropical  climate,  noon. 

B-1 6  Spatial  correlation  of  wind  direction,  coastal  climate,  midnight. 
B-1 7  Spatial  correlation  of  wind  direction,  coastal  climate,  sunrise. 

Bri  8  Spatial  correlation  of  wind  direction,  coastal  climate,  noon. 
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B-3  Spatial  correlation  of  wind  direction,  continental  climate,  noon. 


B-5  Spatial  correlation  of  wind  direction,  arctic  climate,  sunrise. 


B-6  Spatial  correlation  of  wind  direction,  arctic  climate,  noon. 


B-7  Spatial  correlation  of  wind  direction,  desert  climate,  midnight. 


B-8  Spatial  correlation  of  wind  direction,  desert  climate,  sunrise. 
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B-9  Spatial  correlation  of  wind  direction,  desert  climate,  noon. 


B-10  Spatial  correlation  of  wind  direction,  maritime  climate,  midnight. 
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al  correlation  of  wind  direction,  maritime  climate,  sunrise. 
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B-12  Spatial  correlation  of  wind  direction,  maritime  climate,  noon. 


B-13  Spatial  correlation  of  wind  direction,  tropical  climate,  midnight. 
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B-1 4  Spatial  correlation  of  wind  direction,  tropical  climate,  sunrise. 
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B-15  Spatial  correlation  of  wind  direction,  tropical  climate,  noon. 


B-16  Spatial  correlation  of  wind  direction,  coastal  climate,  midnight. 


J-17  Spatial  corralation  of  wind  diraction,  coastal  climat©,  sunrise. 


APPENDIX  C 

SPATIAL  CORRELATION  OF  WIND  SPEED 

C-1  Spatial  correlation  of  wind  speed,  continental  climate,  midnight. 
C-2  Spatial  correlation  of  wind  speed,  continental  climate,  sunrise. 
C-3  Spatial  correlation  of  wind  speed,  continental  climate,  noon. 
C-4  Spatial  correlation  of  wind  speed,  arctic  climate,  midnight. 

C-5  Spatial  correlation  of  wind  speed,  arctic  climate,  sunrise. 

C-6  Spatial  correlation  of  wind  speed,  arctic  climate,  noon. 

C-7  Spatial  correlation  of  wind  speed,  desert  climate,  midnight. 

C-8  Spatial  correlation  of  wind  speed,  desert  climate,  sunrise. 

C-9  Spatial  correlation  of  wind  speed,  desert  climate,  noon. 

C-10  Spatial  correlation  of  wind  speed,  maritime  climate,  midnight. 
C-1 1  Spatial  correlation  of  wind  speed,  maritime  climate,  sunrise. 
C-1 2  Spatial  correlation  of  wind  speed,  maritime  climate,  noon. 

C-1 3  Spatial  correlation  of  wind  speed,  tropical  climate,  midnight. 
C-1 4  Spatial  correlation  of  wind  speed,  tropical  climate,  sunrise. 

C-1 5  Spatial  correlation  of  wind  speed,  tropical  climate,  noon. 

C-1 6  Spatial  correlation  of  wind  speed,  coastal  climate,  midnight. 
C-1 7  Spatial  correlation  of  wind  speed,  coastal  climate,  sunrise. 

C-1 8  Spatial  correlation  of  wind  speed,  coastal  climate,  noon. 
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C-15  Spatial  correlation  of  wind  speed,  tropical  climate,  noon. 
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APPENDIX  D 

TEMPORAL  CORRELATION  OF  TEMPERATURE 

D-1  Temporal  correlation  of  temperature,  continental  climate,  midnight. 
D-2  Temporal  correlation  of  temperature,  continental  climate,  sunrise. 
D-3  Temporal  correlation  of  temperature,  continental  climate,  noon. 
D-4  Temporal  correlation  of  temperature,  arctic  climate,  midnight. 

D-5  Temporal  correlation  of  temperature,  arctic  climate,  sunrise. 

D-6  Temporal  correlation  of  temperature,  arctic  climate,  noon. 

D-7  Temporal  correlation  of  temperature,  desert  climate,  midnight. 

D-8  Temporal  correlation  of  temperature,  desert  climate,  sunrise. 

D-9  Temporal  correlation  of  temperature,  desert  climate,  noon. 

D-10  Temporal  correlation  of  temperature,  maritime  climate,  midnight. 
D-1 1  Temporal  correlation  of  temperature,  maritime  climate,  sunrise. 
D-1 2  Temporal  correlation  of  temperature,  maritime  climate,  noon. 

D-1 3  Temporal  correlation  of  temperature,  tropical  climate,  midnight. 
D-1 4  Temporal  correlation  of  temperature,  tropical  climate,  sunrise. 

D-1 5  Temporal  correlation  of  temperature,  tropical  climate,  noon. 

D-1 6  Temporal  correlation  of  temperature,  coastal  climate,  midnight. 
D-1 7  Temporal  correlation  of  temperature,  coastal  climate,  sunrise. 

D-1 8  Temporal  correlation  of  temperature,  coastal  climate,  noon. 
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D-2  Temporal  correlation  of  temperature,  continental  climate,  sunrise. 
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D-3  Temporal  correlation  of  temperature,  continental  climate,  noon. 
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D-4  Temporal  correlation  of  temperature,  arctic  climate,  midnight. 
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Temporal  correlation  of  temperature,  arctic  climate,  sunrise. 
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D-6  Temporal  correlation  of  temperature,  arctic  climate,  noon. 
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D-7  Temporal  correlation  of  temperature,  desert  climate,  midnight. 
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Temporal  correlation  of  temperature,  desert  climate,  sunrise. 
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D-10  Temporal  correlation  of  temperature,  maritime  climate,  midnight. 
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Temporal  correlation  of  temperature,  maritime  climate,  sunrise. 
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D-12  Temporal  correlation  of  temperature,  maritime  climate,  noon. 
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D-13  Temporal  correlation  of  temperature,  tropical  climate,  midnight. 
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D-14  Temporal  correlation  of  temperature,  tropical  climate,  sunrise. 
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15  Temporal  correlation  of  temperature,  tropical  climate,  noon. 
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D-16  Temporal  correlation  of  temperature,  coastal  climate,  midnight. 


APPENDIX  E 

TEMPORAL  CORRELATION  OF  WIND  DIRECTION 


E-1  Temporal  correlation  of  wind  direction,  continental  climate,  midnight. 
E-2  Temporal  correlation  of  wind  direction,  continental  climate,  sunrise. 
E-3  Temporal  correlation  of  wind  direction,  continental  climate,  noon. 

E-4  Temporal  correlation  of  wind  direction,  arctic  climate,  midnight. 

E-5  Temporal  correlation  of  wind  direction,  arctic  climate,  sunrise. 

E-6  Temporal  correlation  of  wind  direction,  arctic  climate,  noon. 

E-7  Temporal  correlation  of  wind  direction,  desert  climate,  midnight. 

E-8  Temporal  correlation  of  wind  direction,  desert  climate,  sunrise. 

E-9  Temporal  correlation  of  wind  direction,  desert  climate,  noon. 

E-10  Temporal  correlation  of  wind  direction,  maritime  climate,  midnight. 
E-1 1  Temporal  correlation  of  wind  direction,  maritime  climate,  sunrise. 
E-1 2  Temporal  correlation  of  wind  direction,  maritime  climate,  noon. 

E-1 3  Temporal  correlation  of  wind  direction,  tropical  climate,  midnight. 
E-1 4  Temporal  correlation  of  wind  direction,  tropical  climate,  sunrise. 

E-1 5  Temporal  correlation  of  wind  direction,  tropical  climate,  noon. 

E-1 6  Temporal  correlation  of  wind  direction,  coastal  climate,  midnight. 
E-1 7  Temporal  correlation  of  wind  direction,  coastal  climate,  sunrise. 

E-1 8  Temporal  correlation  of  wind  direction,  coastal  climate,  noon. 
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E-4  Temporal  correlation  of  wind  direction,  arctic  climate,  midnight. 
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Temporal  correlation  of  wind  direction,  desert  climate,  midnight. 


E-10  Temporal  correlation  of  wind  direction,  maritime  climate,  midnight. 
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E-12  Temporal  correlation  of  wind  direction,  maritime  climate,  noon. 
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E-13  Temporal  correlation  of  wind  direction,  tropical  climate,  midnight. 


E-14  Temporal  correlation  of  wind  direction,  tropical  climate,  sunrise. 
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E-15  Temporal  correlation  of  wind  direction,  tropical  climate,  noon. 


APPENDIX  F 


TEMPORAL  CORRELATION  OF  WIND  SPEED 

F-1  Temporal  correlation  of  wind  speed,  continental  climate,  midnight. 
F-2  Temporal  correlation  of  wind  speed,  continental  climate,  sunrise. 
F-3  Temporal  correlation  of  wind  speed,  continental  climate,  noon. 

F-4  Temporal  correlation  of  wind  speed,  arctic  climate,  midnight. 

F-5  Temporal  correlation  of  wind  speed,  arctic  climate,  sunrise. 

F-6  Temporal  correlation  of  wind  speed,  arctic  climate,  noon. 

F-7  Temporal  correlation  of  wind  speed,  desert  climate,  midnight. 

F-8  Temporal  correlation  of  wind  speed,  desert  climate,  sunrise. 

F*9  Temporal  correlation  of  wind  speed,  desert  climate,  noon. 

F-10  Temporal  correlation  of  wind  speed,  maritime  climate,  midnight. 
F-1 1  Temporal  correlation  of  wind  speed,  maritime  climate,  sunrise. 
F-1 2  Temporal  correlation  of  wind  speed,  maritime  climate,  noon. 

F-1 3  Temporal  correlation  of  wind  speed,  tropical  climate,  midnight. 
F-1 4  Temporal  correlation  of  wind  speed,  tropical  climate,  sunrise. 

F-1 5  Temporal  correlation  of  wind  speed,  tropical  climate,  noon. 

F-1 6  Temporal  correlation  of  wind  speed,  coastal  climate,  midnight. 
F-1 7  Temporal  correlation  of  wind  speed,  coastal  climate,  sunrise. 

F-1 8  Temporal  correlation  of  wind  speed,  coastal  climate,  noon. 
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5  Temporal  correlation  of  wind  speed,  arctic  climate,  sunrise. 
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F-7  Temporal  correlation  of  wind  speed,  desert  climate,  midnight. 
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F-1 1  Temporal  correlation  of  wind  speed,  maritime  climate,  sunrise. 
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F-12  Temporal  correlation  of  wind  speed,  maritime  climate,  noon. 
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